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Abstract: A method of chaos reduction for Hamiltonian systems is applied to 
control chaotic advection. By adding a small and simple term to the stream 
function of the system, the construction of invariant tori has a stabilization effect 
in the sense that these tori act as barriers to diffusion in phase space and the 
controlled Hamiltonian system exhibits a more regular behaviour. 
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1. INTRODUCTION 

Chaotic advection (or Lagrangian chaos) was in- 
troduced by Aref (1988) to qualify a motion in 
which it is possible to generate chaotic trajecto- 
ries even if the flow is laminar. Such phenomena 
are observed in a wide range of physical systems 
(Ottino (1989)) and has fundamental applications 
for instance in geophysical flow (Behringer, et 
al. (1991); Brown and Smith (1991)). Chaotic 
advection was mostly studied for two dimen- 
sional unsteady flow (Solomon and Gollub (1988); 
Solomon, et al. (2003); Camassa and Wiggins 
(1991)). The advantage of these studies is that it 
uses the theory of dynamical systems and in par- 
ticular the Hamiltonian theory for incompressible 
flows; the phase space being the physical space in 
this case. For these motions, experiments can be 
implemented in a laboratory. 
Even though chaos is sometimes preferable in 



mixing problems (Benzekri, et al.(2006)), we usu- 
ally want to suppress it for instance in chaotic 
advection. Recently, a method to control chaotic 
diffusion in Hamiltonian dynamical systems was 
proposed by (Vittot, et al. (2005); Chandre, et 
al.(2005)). It was shown that it is possible to 
prevent chaotic diffusion by adding a small term 
to the Hamiltonian. This technique was used to 
a model describing anomalous transport in mag- 
netized plasmas, and applied experimentally on a 
beam of electrons produced by a long Travelling 
Wave Tube (Chandre, et al.(2005)). 
Our goal in this paper is to use this method 
of control of Hamiltonian chaos for the prob- 
lem of chaotic advection in a two dimensional 
time periodic flow. We apply this method within 
the framework of an experiment using magneto- 
hydrodynamic technique which shows that parti- 
cle trajectories in a time periodic flow are chaotic. 
More precisely, the experiment consists in an elec- 
tric current passing through a thin layer of an elec- 
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trolytic solution with a free surface. The dynamics 
of passive particles of a time periodic two di- 
mensional and incompressible flow is Hamiltonian. 
The Hamiltonian modeling the flow was derived 
ad hoc in (Solomon and Gollub (1988); Solomon, 
et al. (2003)). A comparison was made with the 
experiment to validate the model. The equation 
of motion of these passive particles comes from 
a Hamiltonian of 1.5 degrees of freedom that we 
study in the limit of weak amplitude oscillations. 
The model proposed to describe this phenomena 
is based on the following streamfunction: 

v E'e(x, y, t) — asin(x + esint) siny, (1) 

where a is the maximal vertical velocity in the 
flow, e is the amplitude of the lateral oscillations 
of the velocity field. The current interacts with 
an alternative magnetic field produced by mag- 
nets below the fluid. A chain of vortices are then 
observed. Time periodic dependence is imposed 
externally with a plunger that oscillates slowly up 
and down, displacing the flow laterally, i.e., in the 
direction perpendicular to the roll axes, giving rise 
to chaotic advection. The same phenomenon is 
observed in Rayleigh-Bcnard convection. For the 
steady state, rolls are formed periodically, ft has 
been shown that by imposing a sinusoidal flow, 
chaotic advection occurs. 

Under the perturbation, i.e. for e ^ 0, the vertical 
heteroclinic connection breaks down and the sta- 
ble and unstable manifold intersect transversely 
thus generating chaotic advection of passive par- 
ticles. 

In this method, we seek for one of the simplest 
perturbations to create barriers around the bro- 
ken separatrix of the integrable case. We choose 
a perturbation depending only on the position 
variables. 

We will use the stream function (1) as the starting 
point to control or reduced the chaotic advection. 

2. LOCAL CONTROL METHOD 

The local control method has been extensively 
described in (Vittot, ct al. (2005)) where the 
corresponding rigorous results were proved. We 
summerize here the results of this paper. For a 
Hamiltonian system with L degrees of freedom, 
the perturbed Hamiltonian is 

H(A,0) = w A + V(A,9), 

where (A, 9) e R L x T L are action-angle like 
variables and u> is a non-resonant vector of M. L . 
Without loss of generality, let us consider a region 
near A = (by translation), since the Hamilto- 
nian is nearly integrable, the perturbation V has 
constant and linear parts in actions of order e, i.e. 

V(A, 9) = ev(0) + ew(0) • A + Q(A, 9), (2) 



where Q is of order 0(\\ A|| 2 ). Note that for e = 0, 
the Hamiltonian H has an invariant torus with 
frequency vector u> at A = for any Q not 
necessarily small. The controlled Hamiltonian is 
then constructed as: 

H c (A,9)=u,-A + V(A,9) + f(9). (3) 

In this situation, the control term / only depends 
on the angle variables and is given by 

f(9) = V(0, 9)-V (-Td V(0, 8), 6) , (4) 

where Y is a linear operator defined MS (\ pseudo- 
inverse of u> ■ dg, i.e. acting on V = J2k Vke ,e k 
as 

YV= Y -^-e Wk . (5) 

Note that / is of order e 2 . This can be seen from 
Eq. (2) since / can be rewritten as 

f(0) = e 2 w(0) • Yd e v - Q {-eYdgv, 6) , 

and Q is quadratic in the actions. For any pertur- 
bation V, Hamiltonian (3) has an invariant torus 
with frequency vector close to u). The equation of 
the torus which is constructed by the control is 

A = -Yd e V(0,6), (6) 

which is of order e since V(0, 9) is of order s. 

In the next section, we will see that the amplitude 
of the control term is small compared with the 
perturbation. 

2.1 Application to the suppression of chaotic 
advection 

In this section, we apply the method previously 
summerized to reduce chaotic transport of passive 
tracers. With this method, we create isolated bar- 
riers of transport. In particular, a barrier created 
between two convection rolls allows one to reduce 
the diffusion across these rolls. Let us first give 
some results in the integrable case. 
For e = 0, the Hamiltonian is integrable, the 
trajectories of advected particles coincide with 
streamlines. The phase space, which is here the 
physical space, is characterized by a chain of rolls 
with separatrices localized at x = mw, with m G 
Z. The fluid is limited by two invariant surfaces 
y = 7T and y = corresponding to the top and 
bottom roll boundaries. 

In the integrable case, as mentionned in (Camassa 
and Wiggins (1991)), the flow is characterized by 
hyperbolic fixed points along the two invariant 
surfaces y = it and y = and localized at x — mir, 
m E Z. These points are joined by a vertical 



heteroclinic connection for which the stable and 
unstable manifolds coincide. 

In order to built a barrier, we select a surface 
which is located around x = 0. However we could 
also create a barrier at x — ir (mod 2ir) 
We map the time-dependent stream function \l/ £ 
given by Eq. (1) into an autonomous Hamilto- 
nian written as H(x,y, E,t) — E + ^e{x, y, r). 
This corresponds to an extension of phase space 
(x,y,E,r), where A = (x, E) and = (y,r) are 
momenta and positions. We assume that at t = 0, 
t(0) = for H . Therefore the equation of motion 
of H are the same as the ones for *I> e since r = t. 
We rewrite the autonomous Hamiltonian in the 
form: 

H(x, y,E,r) = E + asin(a; + esinr) siny, 

= Ho(x,E) + * e (x,y,T), (7) 

and develop the stream function in Taylor series 
around x = 

H(x, y, E,t) = E + ev(y, r) + ew(x, y, t)x 

+ Q(x,E,y, T ), (8) 

where H n (x,E) = E is the integrablc part, 
ev(y,r) = * e (0,y,r) and 
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Fig. 1. Streamlines at (a) t = and (6) t = 
37r/4, and (c) Poincare section of the stream 
function (1). The parameters are a = 0.6 and 
e = 0.63. 

where 

C e (t ) = T sin(e sin t) , 

= E^7^2„+i C0S ( 2 «+ 1 K (11) 

n>0 

The equation of the invariant torus is given by 
Eq. (6) 



ew(x,y,r) =d x ^ € (0,y, r), 

oc 1 

Q(x, E, y,r) = Y, 77-W (°- ^ ^ • 

i=i \ >' 

The frequency vector u) is given by 

ijj = (dH /dx,dHo/dE) = (0, 1) and is resonant. 

In order to compute the operator T given by 
Eq.(5), we rewrite ty e by expanding sin(esinTj) 
and cos(e sin t) as series of Bessel functions of first 
kind: 

* e (x, y,t) = a sin y sin x(Jq + 2 J% n cos 2nf) 

n>l 

+ 2a sin y cos x(S^ J^n+i sin(2« + 

n>l 

where J[ — Ji{e) for I G N, and Ji are Bessel 

functions of the first kind. 

The control term is given by Eq.(4): 

f(y,t) = * e (0,y,t) - * e (-ra y * e ,y,t). 

(9) 

Since ^ e does not depend on E, we only have to 
compute Td y ^' e (0 7 y, t) and then 

f(y, t) = a sin(e sin t) siny 

—a sin [— acosyC e (t) + esint] sin?|10) 



x = a cos yC e (t) . (12) 
The controlled stream function is given by: 

^ c (x, y, t) — a siny (sin (x + esint) + sin(esint) 
— sin [— acosyC e (t) + esinf]) . (13) 

We notice that with a small modification of the 
stream function, provide that a and e are such 
that | ae \< 1 , there is an exact formula giving 
the equation of the torus. This invariant torus 
suppress the chaotic transport from roll to roll 
and then along the channel. 

To illustrate numerically these results, we consider 
fixed values of the parameter a and e. Similar 
results hold for other values of the parameters a 
and e. 

Streamlines of the streamfunction (1) are depicted 
on Fig. 1(a) and (6) at two different times t = 
and t — 37r/4 respectively and for e = 0.63 and 
a = 0.6. We observe closed curves corresponding 
to lateral oscillations in the x-direction with a pe- 
riodic displacement of — e sin t. For the same values 
of e and a, a Poincare section of the dynamics of 
the streamfunction (1) for initial conditions with- 
out control are represented in Fig. 1(c). It shows 
that the transport along the channel is greatly 
enhanced. The phase space is characterized by 
regular (quasiperiodic) trajectories and a chaotic 
region around and between the rolls. 
When we add the control term, the streamlines 




Fig. 2. Streamlines at (a) t = and (6) t = 
37r/4, and (c) Poincare section of the stream 
function (13). The parameters are a = 0.6 
and e = 0.63. 

of the controlled stream function (13) are then 
slightly modified (non-uniformly in y) as it can 
be seen on Fig. 2(a) and (6) for e = 0.63 and 
a = 0.6 at two different times t = and t = 37r/4 
respectively. Moreover, the displacement of these 
rolls remains parallel to the x-axis as it is for 
the stream function given by Eq. (1). A Poincare 
section of the dynamics of passive tracers in a 
flow described by the stream function \t c given 
by Eq. (13) is represented in Fig. 2(c). We notice 
first that there are invariant surfaces which have 
been created around x — (mod 2tt) (bold curves) 
which prevent the diffusion of passive particles 
along the channel. Moreover we observe a regular- 
isation of the dynamics whithin the cell bounded 
by the two barriers x — and it (mod it). 

In Fig. 3, we depict a numerical simulation of the 
dynamics of a dye in the fluid. The left column 
shows the evolution of the tracers for the stream 
function \P e given by Eq. (1). The right column 
shows that the control term regularizes the dy- 
namics and prevents any spread of a dye within 
a cell which is limited by two barriers created by 
the stream function ^f c given by Eq. (13). 

For | ae \« 1, the control term (13) can be 
simplified such as: 



f s {y,t) = -— sin2ycos(esint)C £ (t). (14) 



We remark that the control term f s is of order a 2 e 
and does not depend on x as it is expected from 
the method. 

In order to test the robustness of the method and 
to try an experimentally more tractable pertur- 
bation, we truncate the series given by Eq. (11) 
by considering only the first term of the series 
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Fig. 3. Numerical simulation of the dynamics 
of a dye at t = 30, t = 50, t = 70 and 
t = 140 (from top to bottom) : left column 
for the stream function (1) and right column 
for the stream function (13). The parameters 
are a = 0.6 and e = 0.63. 

C e (t) in Eq. (11) and cos(esini) which leads to 
the following stream function 

\E' s (x, y,t) — a sin(x + e sin t) siny 

-a 2 Jiaia 2y cost. (15) 

The perturbation has a norm sup ]t | V(x,y,t) | = 
e and the norm of the the simplified control term 
fs is 



sup | f s (y,t) \=a z Jl 
y,t 



(16) 



Therefore the ratio r between the norm of the 
control term and the perturbation is given by 
r w ae. 

The dynamics of the stream function ^ s is de- 
picted in Fig. 4 for e = 0.63 and a = 0.6. We 
see that invariant surface has been created around 
x = (mod 2tt) and that the control term still re- 
duces significantly the chaotic advection but there 
are some advected particles along the channel. We 
observe that an other invariant surface has been 
created around x — ir. This is due to the fact 
that the flow given by the stream function \P S is 
invariant under the symmetry 

t — > t, x — > x + it, y — *■ y + it. 

It results from this symmetry, a regularisation of 
the dynamics within the cell bounded by the two 
barriers x = and it (mod 2tt). This symmetry is 
an approximate one ( up to order | a 2 e |) for the 
stream function ty c . Thus the control term is able 
to regularize its dynamics also in this region. 



Fig. 4. Poincare section of the stream function 
(15) . The parameters are a — 0.6 and e = 
0.63. 

In order to see more clearly the effect of the 
control term, we study the diffusion properties 
of the system. The mean square displacement 
< r 2 (t) > of a distribution of Ai particles (of 
order 3000) is computed as a function of time: 

M 

<r 2 (t)>=—Y J W^-^)\\\ (17) 

i=i 

where Xj(f), i = 1, M is the position of the i-th 
particle at time t obtained by integrating Hamil- 
ton's equations with initial conditions Xi(0). We 
find that < r 2 (t) > grows linearly for the consid- 
ered time interval (see Fig. 5(a)). The transport is 
then assumed to be described as normal diffusion 
and the corresponding diffusion coefficient can be 
determined from the slope of < r 2 (t) > versus t: 

D=lim <r2(f)> . (18) 

t^QO t 

Figure 5(6) shows the values of D as a function of 
e with and without control term (15) determined 
from the mean square displacement for t > 1000. 
We remark that the diffusion coefficient with the 
control term (15) is significantly smaller than the 
uncontrolled case. 

3. CONCLUSION 

We presented in this paper an application of a 
recently developed control technique to Rayleigh- 
Benard convection. We derived analytically the 
expansion of the control term required to reduce 
chaotic advection present in the original prob- 
lem. Using Poincare sections, and tracking the 
dynamics of a dye, we showed the efficiency of 
the method. The slightly deformed streamlines 
due to the control terms prevent any large spread 
through the channel but confine the material lines 
inside some deformed rolls. The chaotic behavior 
of the flow was significantly reduced by taking 
only the first order of control term. 
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Fig. 5. (a) Mean square dispacement < r 2 (t) > 
versus time t obtained for Hamiltonian 
(1) with three different values of e = 
0.12,0.18,0.2 and (6) Diffusion coefficient D 
versus e obtained for stream function (1) (cir- 
cle) and stream function (15) (square) 
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